Extraordinary magnetoresistance in graphite: experimental evidence for the 

time-reversal symmetry breaking. 
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We report an unusual highly-anisotropic in-plane magneto-resistance (MR) in graphite. In a 
certain current direction MR is negative and linear in low fields with a crossover to a positive MR at 
higher fields, while in a perpendicular current direction we observe a giant super-linear positive MR. 
These extraordinary MRs are respectively explained by a hopping magneto-conductance via non-zero 
angular momentum orbit als, and by the magneto-conductance of inhomogeneous media. The linear 
negative orbital MR is a unique signature of the broken time-reversal symmetry (TRS) in graphite. 
While some local paramagnetic centers could be responsible for the broken TRS, the observed large 
diamagnetism suggests a more intriguing mechanism of this breaking, involving superconducting 
clusters with unconventional (chiral) order parameters and spontaneously generated normal-state 
current loops in graphite. 
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The transverse magneto-resistance (MR) in a mag- 
netic field B along the z-axis defined as MR = 
o'xxO'xx{^)/{o'xx + ^xy) " 1 tcrms of the conductivity 
tensor (Jik{B)^ is a sensitive probe of the electron trans- 
port mechanism in doped semiconductors and metals. In 
an isotropic medium with the Bloch electrons it is posi- 
tive and quadratic in B, except the open Fermi surfaces, 
where for some specific directions of B the positive MR 
is linear in B. In the hopping regime with localized carri- 
ers MR is caused by a strong magnetic field dependence 
of the exponential asymptotic of bound state wave func- 
tions at a remote distance from a donor (or an acceptor), 
which leads to a significant positive MR, which is also 
quadratic in low magnetic fields [2^. 

However, there is a substantial number of semiconduc- 
tors and semimetals (e. g. bismuth) where open Fermi 
surfaces are not feasible, but the positive MR is linear. 
One of the theoretical possibilities for such a phenomenon 
is a so-called quantum magnetoresistance in semimetals 
having pockets of the Fermi surface with a small or even 
zero effective mass (as the Dirac fermions in graphite and 
graphene), which might be in the ultra-quantum limit at 
rather low magnetic fields [3]. 

On the other hand, there is anomalous negative MR 
(NMR) first observed in some hopping systems, for in- 
stance in amorphous germanium and silicon. More re- 
cently a linear NMR has been observed in the longitudi- 
nal c-axis interlayer current in graphite at high magnetic 
fields assigned to a growing population of the zero-energy 
Landau level of quasi-two-dimensional Dirac fermions 
with the increasing magnetic field [4]. 

Here we report a low- field transverse MR in bulk quasi- 
two-dimensional graphite samples, which is giant ( 1400% 
at 2K and B = 200 mT) and positive in one in-plane cur- 
rent direction, and linear and negative (below 80 mT) 
in the other direction, when the magnetic field is ap- 
plied perpendicular to the planes. These observations 
are quantitatively explained in the framework of the mag- 
netotransport of inhomogeneous low-carrier semiconduc- 



tors [5] and of a hopping magneto-conductance via non- 
zero angular momentum orbitals [6]. 

We have measured 10 commercially available HOPG 
samples with different mosaicity, characterized by 
FWHM of X-ray rocking curves, and the room temper- 
ature out-of-plane/basal plane zero-field resistivity ratio 
Pc/ pah ranging from 10^ to 10^. X-ray diffraction mea- 
surements revealed a characteristic hexagonal graphite 
structure in the Bernal [ABAB) stacking configuration, 
with no signature for other phases. The obtained crys- 
tal lattice parameters are a = 2.48A and c = 6.7lA. 
The results given in this letter were obtained on most 
anisotropic {pc/ Ph = 2 x 10^) HOPG samples from Union 
Carbide Co. and here shown for the sample with = 
S/il^cm, pc = ir^cm, and FWHM = 0.7^ (the sample size 
i^lxwxt = 2.5 X 2.5 x 0.5mm^). The magnetic field was 
applied parallel to the hexagonal c-axis {B \\ c \\ t), and 
Pah{B^T) were recorded using the van der Pauw geome- 
try (see insets in Fig. 1) sweeping the field between a tar- 
get negative and positive values using Quantum Design 
and Janis 9 T-magnet He-4 cryostats. For the measure- 
ments, silver past electrodes were placed on the sample 
surface, while the resistivity values were obtained in a ge- 
ometry with a uniform current distribution through the 
sample cross section. All resistance measurements were 
performed in the Ohmic regime. Complementary mag- 
netization measurements M{B^T) were carried out by 
means of a SQUID (Quantum Design) magnetometer. 

The graphite samples described above show qualita- 
tively different temperature dependence of the in-plane 
zero-field resistance for the current in two perpendicu- 
lar in-plane directions, Fig.Q. In one direction (called 
here Y) the resistance is metallic-like while in the other 
direction (X) it is insulating-like in a wide temperature 
interval above 2K. Since a sufficiently strong magnetic 
field makes graphene planes electrically isotropic, Fig.([T]), 
the anisotropy is of electronic origin. Most probably it 
is associated with the inhomogeneous carrier-density dis- 
tribution, such that metallic clusters are partially over- 
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FIG. 1: (Color online) Upper panel: Insulating- like (X) and 
metallic- like (Y) zero-field in-plane resistance of graphite in 
two perpendicular current directions. The resistance becomes 
isotropic and insulating-like in high magnetic fields. Lower 
panel: well doped metallic clusters (dark regions) overlapped 
in Y-direction in an insulating matrix (light regions). 



lapped in Y-direction while they are separated by poorly 
doped insulating regions in X-direction, as illustrated 
in the lower panel of Fig.([T]). In high magnetic fields 
about 1 Tesla and higher resistance becomes isotropic 
and insulating-like in all directions, Fig.([T]), because of a 
giant positive MR of the metallic clusters and suppres- 
sion of possible superconductivity, as discussed below. 

When a relatively weak magnetic field is applied per- 
pendicular to the planes, the MR along the metallic Y 
direction appears huge and positive, Fig.(|2]). A giant pos- 
itive MR is naturally expected in doped graphite with 
virtually massless Dirac fermions 0] since the parame- 
ter j3 = uJcT becomes large already in the mT-region of 
the field (here uOc and r are the Larmor frequency and 
the scattering time, respectively). However, it is neither 
quadratic as in the Boltzmann theory, nor linear in B as 
in the quantum magnetoresistance [3]. We suggest that 
inhomogeneities are responsible for the strong departure 
from these regimes. They lead to a radical rearrangement 
of the current flow pattern changing the magnetic field 
dependence of the transverse conductivity in the strong 
magnetic field, {3 > 1 [5]. Importantly, when cjxy ^ (Jyy^ 
then even relatively small inhomogeneities in the carrier 
density lead to the MR proportional to B^^^ [5[. In fact, 
MRy = {B / BinY^^ fits nicely the observed magnetic 
field dependence of the magnetoresistance in the metallic 



FIG. 2: (Color online) Magnetoresistance in the metallic (Y) 
direction at different temperatures (symbols). 



current direction, Fig.(|3]) with a single scaling parameter 
Bin depending on the fluctuations in the carrier density 

The MR in the insulating direction is even more 
anomalous. Fig.®. It is linear and negative in very low 
fields and positive and superlinear in higher fields above 
the crossover point. The insulating-like temperature de- 
pendence of the zero-field resistance, Fig.([T]), in this di- 
rection supports the view that the metallic clusters are 
virtually nonoverlaped along X, so that the resistance is 
the sum of the hopping resistance, Rh of the insulating 
regions with low doping, and the metallic resistance Km- 
Hence, the magnetoresistance in X direction is expressed 
via metallic MRy and the hopping magnetoresistance, 
MRh of insulating layers as 



MRx 



MRy MRh 



(1) 



where r = Rh{0) / Rm{0) is the ratio of the zero-field in- 
sulating resistance to the zero-field metallic resistance. 
This ratio is about r ^ 8.5 at 2K, Fig. ([I]). 

In inorganic and organic insulators lattice defects such 
as vacancies, inter st it ials, excess atoms or ions and other 
"impurities" often localize carriers with a finite momen- 
tum rather than in the zero- momentum s-states [8]. To 
model the hopping MRh we apply a recent development 
of the conventional theory of hopping magnetoresistance 
for hopping via non-zero orbital momentum states [6|. 
Quite remarkably this renders a large weak-field hopping 
MRh, which is negative. The theory [6] has accounted 
for the unconventional MRs in some organic insulators (a 
so-called OMAR ) observed in about the same range of 
the magnetic field [9] as in our graphite samples. More- 
over, if the orbital degeneracy is lifted due to a broken 
time-reversal symmetry with or without net magnetiza- 
tion, the negative MRh is linear in 5 [6], 
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FIG. 3: (Color online) Upper panel: Magnetoresistance in the 
insulating (X) direction at different temperatures (symbols). 
Lower panel shows the low- field negative MR fitted as MRx = 
-cB with c 0.0029 1/mT (T=2K), c 0.00115 1/mT 
(T=20K), c = 0.00053 1/mT (T=30K), c = 0.00028 1/mT 
(T=40K), indicating the time-reversal symmetry breaking. 



Here m is the magnetic quantum number of the localized 
state, is the inverse zero-field localization length, r and 
p are an average hopping range and its projection on the 
plane perpendicular to the magnetic field, respectively, 
and b = B / Bq is the reduced magnetic field with Bq = 

For the s-wave bound states with m = Eq.([2]) is 
the textbook result [2] predicting the positive hopping 
MRh. On the contrary, for hopping via the orbitals 
with nonzero orbital momentum MRh is negative and 
linear in relatively small B, if m = 1,2,3,.... The linear 
NMR is caused by the expansion of the localized states 
with the positive m due to a linear magnetic lowering of 
their ionization energy by hjjJcm/2. It dominates in the 
wide region of realistic impurity densities for any nonzero 
m because of the small numerical factor (1/48) in the 
quadratic term in the exponent of Eq.([2]). 

Using MRh ^ ex.p{—B/Bh) — 1 in Eq.([T]) fits remark- 
ably well the low-temperature MR in the insulating cur- 
rent direction with a single scaling parameter Bh pro- 
portional to the localized-state ionization energy, Fig.(|3]). 
The same Eq.([T]) describes reasonably well the negative 
to positive MR crossover even at temperatures as high as 
50K. With increasing temperature deeper localized states 
with a higher ionization energy become accessible for the 
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FIG. 4: (Color online) Magnetoresistance in the insulating 
(X) direction at T = 2K (symbols) fitted by Eq.p with r = 
8.5 and Bh — 110 mT (line). Inset: magnetoresistance in 
the metallic (Y) direction at T = 2K (symbols) fitted by 
(B/Bin)^^^ with Bin = 26.9 mT (line). 



hopping conductance, so that the slope c = l/Bh of the 
linear negative MR drops, as observed. Fig. ([3]). 

While the outhned model with one or two scaling pa- 
rameters provides accurate agreement with our experi- 
mental observations, there are other theoretical mecha- 
nisms of NMR. Importantly, there is no NMR and vir- 
tually no MR in the same range of the magnetic field 
parallel to the planes, which rules out a spin origin of the 
observed MRs. Also there are no quantum magnetic os- 
cillations at high temperatures, where unusual MRs are 
still observed, Figs. (j2|3p . so they are not related to the 
Landau quantization. The weak localization gives NMR 
which is often almost linear in a certain field range. But 
such NMR smoothly evolves from a sub-linear magnetic 
field dependence at lower temperatures to super-linear 
field dependence at higher temperatures, while we ob- 
serve a perfectly linear NMR at low 5 in a wide temper- 
ature range, Fig.(j3j). A strong NMR exists in the classical 
two-dimensional electron gas due to freely circling elec- 
trons, which are not taken into account by the Boltzmann 
approach. However, it is parabolic rather than linear 
at low fields [lo|- The parabolic orbital NMR has been 
also predicted by the gauge theory in two-dimensional 
strongly-correlated doped Mott insulators [11]. 

On the other hand, if there is no time-reversal sym- 
metry breaking, the states with the opposite direction of 
the orbital angular momentum, m and — m, are degener- 
ate, so that the linear term in the MRh, Eq.([2]) cancels. 
But ions that carry a magnetic moment break the time- 
reversal symmetry and split m and — m states. Such 
zero-field splitting (ZFS) gives preference to the hopping 
via orbitals with a lower ionization energy (positive m) 
providing the negative linear MR [6]. In fact, we mea- 
sured ferromagnetic M{B) hysteresis loop for B paral- 
lel to graphene planes in the electrically homogeneous 
graphite (not shown here, see Ref. [l2|) demonstrating a 
metallic-like zero-field resistance in all directions at tem- 
peratures below 50 K. MR of those samples can be well 



4 



1 ' 1 ' 1 ' 1 ' 

^'^fev 0.5 
i 00 






B II planes 
T = 2 K 


bj, B (mT) 



-600 -400 -200 200 400 600 

B (mT) 



FIG. 5: Diamagnetic magnetization of electrically inhomo- 
geneous graphite with NMR. Inset shows the magnetization 
after subtracting the initial (B < Bk) linear diamagnetic con- 
tribution, Md{B) = xB with X = -3.4 x 10~^ emu/g mT. 



fitted by B"^^^ law as in Fig.(j4]). 

However, we do not observe a global ferromagnetism 
in the bulk electrically inhomogeneous samples with the 
linear NMR, Fig. ([3]). In a sharp contrast with the elec- 
trically homogeneous samples the electrically inhomoge- 
neous samples show a large diamagnetic response almost 
linear in Fig.(j5j). Naturally, some local paramag- 
netic centers responsible for the broken TRS could be 
found also in the inhomogeneous graphite, with their 
magnetic response overwhelmed by the large diamag- 
netism. However the observed large diamagnetism itself 
suggests a more intriguing mechanism of TRS breaking 
in our inhomogeneous samples, such as superconducting 
clusters 13 1 with an unconventional (chiral) order pa- 
rameter 14,[l5|. There is a kink in the field dependence 



of the diamagnetic magnetization of the inhomogeneous 
samples at B^ ~ 0.2 T, Fig. ([5]), resembling the behav- 
ior of type-II superconductors in magnetic fields exceed- 
ing the lower critical field Bd- Supporting this possi- 
bility the electrically inhomogeneous samples are becom- 
ing perfectly homogeneous insulators in sufficiently high 
magnetic fields, which could suppress the superconduc- 
tivity, Fig.([T]). Also in some 2D lattices there is the bro- 
ken time-reversal symmetry (i.e. some internal magnetic 
ordering) in the normal state with a spontaneous quan- 
tum Hall effect but without any net magnetic flux at high 
temperatures [16]. At lower temperatures a pair of spon- 
taneously generated current loops in adjacent graphene 
layers, having odd-parity with respect to the two layers, 
could also break TRS |17| . 

In conclusion, we presented the highly-anisotropic in- 
plane magnetoresistance in graphite, which, in some cur- 
rent direction, is negative and linear in low fields with 
the crossover to the positive MR at higher fields, while in 
the perpendicular current direction we observed the giant 
super-linear positive MR. The extraordinary MRs have 
been explained by the hopping magneto-conductance 
via non-zero angular momentum orbit als, and by the 
magneto-conductance of inhomogeneous media. More 
generally our findings point to an inhomogeneous doping 
and a semiconducting gap in graphite (see also Ref.[18]). 
The linear orbital NMR is a unique signature of the time- 
reversal symmetry breaking. Combined with the large 
diamagnetism it points towards superconducting clus- 
ters with unconventional order parameter and normal 
state spontaneous currents in electrically inhomogeneous 
graphite samples. 
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